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A search for the process e+e“ —>■ r]' has been performed with the SND detector 
at the VEPP-2000 e'’'e“ collider. The data were accumulated at the center-of-mass 
energy of 957.78 ± 0.06 MeV with an integrated luminosity of about 2.9 pb“^. For 
reconstruction of the rj' meson five decay chains have been used: r]' —>■ r]TT~^TT~ 
followed by the r] decays to 77 and Svr^, and rj' r/vr^vr^ followed by the 7 decays to 
TT+TT^vr^, 77 , and Svr*^. As a result, the upper limit has been set on the rj' electronic 
width: r^/^e+e- < 0.0020 eV at the 90% confidence level. 

PACS numbers: 13.20.Jf, 13.40.Gp, 13.66.Bc, 14.40.Be 


I. INTRODUCTION 

This article is devoted to a search for the rare leptonic decay r]' —>■ In the Stan¬ 

dard model this decay proceeds through the two-photon intermediate state (Fig. [T]) and, 
therefore, is suppressed as o? relative to the r! two-photon decay. An additional suppression 
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FIG. 1: The diagram for the rj' —?> e"*“e decay. 


of {me/m'^Y arises from helicity conservation. The imaginary part of the p' e'^e~ decay 
amplitnde can be expressed in terms of the known two-photon width r{'r]' —)■ 77 ). Neglecting 
the real part of the amplitnde one can obtain the model-independent lower limit (nnitary 


limit) on the decay probability B{r]' e+e ) > 3.8 x 10 [ll. Calcnlation of the real part 
reqnires knowledge of the transition form factor F^ql^qY) for the 7 * 7 * —)■ r]' vertex, where 
4 a„d ,1 are the photoa vatuahties the loop. The real part may merease the decay 
probability by a factor of 3-5 as compared with the nnitary limit [2|, 13[ . Due to the small 
probability, the e’''e“ decay may be sensitive to contributions not described by the 

Standard Model 5|. 

The strictest limit on the decay branching fraction B{ri' e~^e~) < 1.2 x 10“® joj at 
the 90% confidence level (CL) was set recently in the experiments with the CMD-3 detector 
at the VEPP-2000 e^e“ collider jsl. In this experiment the technique of using the inverse 
reaction e+e“ —?• rj' for a measurement of B{ri' —)■ e+e“) proposed in Ref. was applied. 
The cross section of the —?• 7 ' reaction at the center-of-mass (c.m.) energy E = rUrf'C^ 

is equal to 


(To = ^-5(7' -)■ e+e ). ( 1 ) 

In this paper we present the results of the search for the rj' e^e~ decay in the ex¬ 
periments with the SND detector at the VEPP-2000 collider. The SND data used in this 
analysis were collected simultaneously with the CMD-3 data mentioned above. 
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II. DETECTOR AND EXPERIMENT 


The SND detector is described in detail elsewhere {q]. This is a nonmagnetic detector, the 
main part of which is a three-layer spherical electromagnetic calorimeter based on Nal(Tl) 
crystals. The solid angle covered by the calorimeter is 90% of dvr. Its energy resolution 
for photons is (Te/E = 4.2%/ ^E{GeY), and the angular resolution is about 1.5°. The 
directions of charged particles are measured by a tracking system, which consists of the 
9-layer drift chamber and the proportional chamber with readout from cathode strips. The 
tracking system covers a solid angle of 94% of 47 r. The calorimeter is surrounded by a muon 
system, which is used, in particular, for cosmic-background suppression. 

Data used in this analysis with an integrated luminosity of about 2.9 pb~^ were accu¬ 


mulated in 2013 at the c.m. energy close to m„/c^ = 957.78 ± 0.06 MeV 


10|. During the 


data taking period the beam energy was monitored with an absolute accuracy of about 60 
keV by the Back-scattering-laser-light system ^ . The data taking conditions are described 
in detail in Ref. 6 (]. The average value of the c.m. energy is E^^a = 957.68 ± 0.060 MeV; 
its spread is cte^^ = 0.246 ± 0.030 MeV. To obtain the cross section of rj' production in 


the real experimental conditions we have to take into account t 
the initial state, and the energy spread. This was done in Ref. 


le radiative corrections to 
6 j. As the collider energy 


spread (FWHM = 0.590 MeV) is significantly larger than the rj' width T.,,/ = (0.198 ± 0.009) 
MeV lO], the resulting cross section is proportional to the electronic width 


a™(nb) = (6.38±0.23)r,,^e+e-(eV). 


( 2 ) 


It should be noted that the radiative corrections and the energy spread lead to a reduction 
of the cross section compared to the Born one (Eq. ([I])) by a factor of four. 

The search for the process e’''e“ —)■ rj' is performed in Eve decay chains: rj' —)■ ri7i~^7r~ 
with the r] decays to 77 and 37 r°, and t]' — rjTr^TT^ with the rj decays to 7 r+ 7 r“ 7 r°, 77 and 37 r°. 
For luminosity normalization of events with charged particles in the final state the large- 
angle Bhabha scattering is used, while for events containing only photons the luminosity 
is measured using the two-photon annihilation e+e“ 77 . The corresponding integrated 
luminosities are measured to be Lee = 2.91 pb“^ and = 2.82 pb“^. The difference 
between these values, about 3%, gives us a conservative estimate of the luminosity systematic 
uncertainty. 
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FIG. 2: The spectrum of the normalized total energy deposition in the calorimeter for simulated 
events of the process e'’'e“ r]' ^ 7r'^Tr~r], rj —)■ 77 . The arrows indicate the boundaries of the 

selection cut 0.55 < E'cai/.^'cm < 0.9. 

Using different normalizations allows to partly cancel systematic uncertainties associ¬ 
ated with hardware event selection, charged track reconstruction, and beam-generated extra 
tracks. 


III. EVENT SELECTION 
A. Decay chain rj' vr+Tr^ry, rj 77 

This f]' decay channel has the largest probability, about 17%, and the lowest multiplicity 
among the channels studied in this work. Because of the small multiplicity the background 
for this channel arises from almost all e’''e“ annihilation processes. In background processes 
with a small number of photons, as e~^e~{'y) or —)■ 7r’''7r“(7), additional fake 

photons appear as a result of splitting of electromagnetic showers, nuclear interaction of 
pions in the calorimeter, or superimposing beam-generated background. 

At the hrst stage events with two charged particles originating from the interaction re¬ 
gion and two photons are selected. The muon system veto is applied to reject cosmic-ray 
background. The charged particle tracks are htted into a common vertex. Their polar angles 
must be in the range 40° < 6 < 140°. To suppress background from collinear two-body pro¬ 
cesses, mainly from e~^e~, the azimuthal angles of the charged particles are required 
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FIG. 3: Two-dimensional distribution of parameter xfj over Schar, 2+3 for data events (top) and 
simulated events of the e+e“ r]' ^ tt+tt”?/, rj ^ process (bottom). The rectangle in the 
bottom left corner of the plot corresponds to the selection criteria < 30 and T^char, 2+3 < 60 
MeV. 

to satisfy the condition |180° — |0i — 02|| > 10°. The background from nj events 

with photon conversion into an e^e“-pair is rejected by the condition 'ipcc > 20°, where "^cc is 
the open angle between the charged particles. To remove events with fake photons from pion 
nuclear interactions in the calorimeter, the condition on the minimal open angle between 
a charged particle and photon 'ipc-y > 20° is applied. The specific feature of this channel 
is a large energy deposition in the calorimeter T^cai- The distribution of this parameter for 
simulated events of the process e’''e“ —)■ r/' —)■ 7i~^7r~ri —)■ 7 r^ 7 r “27 is shown in Fig. |2l The 
arrows indicate the boundaries of the condition used 0.55 < Ecai/Ecm < 0.9. 

For events passing preliminary selection the kinematic £t to the e’''e“ —)■ hypoth¬ 

esis is performed. The quality of the fit is characterized by the parameter Another 
important parameter used for the final selection is the sum of energy depositions of charged 
particles in the second and third layers of the calorimeter i? 2 + 3 ,char- Since pions in the pro¬ 
cess under study are soft, they stop predominantly in the first layer of the calorimeter. The 
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two-dimensional distributions of the parameters and -E 2 -i- 3 ,char for data events and simu¬ 
lated events of the process under study are shown at Fig. [3l The rectangle in the bottom 
left corner corresponds to the selection criteria applied: < 30 and ii^char, 2+3 < 60 MeV. 

No data events are selected with the selection criteria described above. The detection 
efficiency for —)■ p' —)■ rj ^ 2^ events is determined using Monte Carlo (MC) 

simulation to be (12.2 ± 1.2)%. The quoted error is systematic. For its estimation we 
use the results of the study of data-MC simulation difference in the measurement of the 

I 2 I at Ecm > 1.2 GeV. 


e’^e —>■ TT^TT T] cross section 


The dominant sources of background for this decay mode after applying the selection 
criteria are the processes e’''e“ —)■ 177 , 7 —)■ 7 r+ 7 r“ 7 r° and e^e“ Tr+Tr^vr^Tr^. The number of 
background events is estimated using MC simulation to be 0.7±0.1 and 0.10±0.05 for the first 
and second processes, respectively. The values of the e+e“ rj'y and 7 r+ 7 r“ 7 r° 7 r° 


cross sections were taken from the measurements 


13l-ll5|. The background can be also 


estimated from the two-dimensional distribution shown in Fig. fusing the assumption that 
the xfj and F^ 2 + 3 ,char distribution are independent. The number of background events in the 
signal rectangle < 30, F^ 2 + 3 ,char < 60 MeV) is estimated as n 2 U 3 /n 4 ^ 1 ± 1, where n 2 , 
Us and are the numbers of data events in the regions (30 < < 60, i? 2 + 3 ,char < 60 MeV), 

(x^ < 30, E 2 + 3 ,char > 60 MeV) and (30 < x^ < 60, F; 2 + 3 ,char > 60 MeV), respectively. 

There is also the nonresonant reaction e+e“ tt+tt”// having the same hnal state as 
the process under study. This reaction proceeds through the pr] intermediate state and, 
therefore, is suppressed due to the small phase space of the final particles. Interpola ting the 
result of the £t to the e’''e“ cross section measured at higher energies 12 |], we 

estimate that the nonresonant cross section at E^ra = 960 MeV is about 1.7 pb. Assuming 
the same detection efficiencies for the resonant and nonresonant processes, the nonresonant 
contribution is estimated to be 0.2 events. 


B. Decay chain t}' tt+tt rj, t] ^ Stt® 


In this decay mode the following selection criteria are used. An event must contain 
two charged particles originating from the interaction region and six photons. We require 
the muon-system veto, |180 — |0i — 02 || > 10°, 'i/’c 7 > 20°, 0.5 < F^cai/-Fcm < 0.9, and 

-^2+3,char 


< 90 MeV. 
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FIG. 4: The two-dimensional distribution of the parameters and 473,^0 for data events (top) 
and simulated rj' —>■ vr+vr^Ty, ?/ —)> Svr^ events (bottom). The rectangle corresponds to the selection 
criteria used: < 50 and 500 < < 600 MeV/c^. 

For selected events the kinematic fit is performed to the hypothesis e+e“ —)■ 7 r+ 7 r“ 37 r°. 
The two-dimensional distributions of of the kinematic £t (Xg^o) versus the three 
invariant mass (Mg^o) for data events and simulated events of e^e“ —)■ ry' —)■ 7 r+ 7 r“ry, rj —)■ Svr^ 
process are shown in Fig. IH The following cuts on these parameters are used: Xg^o < 50 
and 500 < Mg^rO < 600 MeV/c^. No data events satisfying the selection criteria applied are 
found. The detection efficiency for e+e“ —)■ ry' ^ 7i~^7r~ri, ry —)■ Svr® events determined using 
MC simulation is (7.5 ± 0.8)%. The quoted error is estimated according to Ref. jl 2 |. 

It is necessary to note that the same hnal state can be obtained in the other 

decay chain rj' —>■ n^7r^T],T] —>■ The distribution of the three- 7 r° invariant mass for 

this decay channel is shown in Fig. [5l 

It is seen that a significant part of t]' —)■ Tr^vr'^ry, ry —)■ Tr+Tr^vr® events satishes the condition 
500 < Mg^rO < 600 MeV/c^. The detection efficiency for this channel calculated using MC 
simulation is (4.9±0.5)%. We could not increase the detection efficiency for rj' —)■ n^TT^r], rj —)■ 
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FIG. 5: The distribution of the three-invariant mass for the decay channels: rj' —>■ 7r^7r“r/,r/ —>■ 
(narrow distribution) and rj' —)> r/ (wide distribution), obtained using MC 

simulation. 

7r'''7r“7r° events by using conditions on parameters specific for this decay mode, for example, 
the TT+TT^TT^ invariant mass instead of Ma^o. 

The dominant background source for the Tr+Tr^Tr^Tr^Tr^ hnal state is the process e+e“ —)■ 
7r+7r“7r°7r°. Additional fake photons can appear as a result of nuclear interaction of charged 
pions or beam background. The number of background events obtained using MC simulation 
is 2.7 ± 0.5. Since appearance of two fake photons is needed, the simulation can be used 
only for rough estimation of the background level. The background can be also estimated 
using the distribution for data. Based on four observed data events with < 50 in 
Fig. [Hand assuming a linear background distribution, we estimate the background in 
the interval 500 < < 600 MeV/c^ to be 2 ± 1 events. The nonresonant background 

from the e^e“ process discussed above in Sec. IIII Al is about 0.1 events in this 

decay mode. 


C. Decay chain rj' —)• rj 77 


The event selection in this decay mode is performed in two steps. At the hrst stage 
six-photon events containing no tracks in the drift chamber are selected. Each photon is 
required to have the transverse energy distribution in the calorimeter consistent with the 
distribution for an electromagnetic shower 16[. The total energy deposition Ecai and the 
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FIG. 6 : The distribution for data events (solid histogram) and simulated e+e rj' ^ 

2TT^r] —>■ 67 events (dashed histogram). 


event momentum P^ai calculated using energy depositions in the calorimeter crystals must 
satisfy the following conditions: 

0.7 < Ecal/Ecra < 1-2, cPca.\/Ecm < 0.3, Pcal/-Ecm “ cPca.l/E^m > 0.7. (3) 


To reject cosmic-ray background the muon-system veto is required. 

For events passing initial selection the kinematic £t to the e’''e“ —>■ —)■ 67 

hypothesis is performed. The quality of the £t is characterized by the parameter x^-n-OirO- 
The distributions of this parameter for data events and simulated events of the e’''e“ —>■ 
rj' —)■ Tr^TT^T] —)■ 67 process are shown in Fig. El The condition X^irO-n-o < 15 is applied. 

No data events satisfying the criteria described above have been found. The detection 
efficiency for the e’''e“ —)■ —)■ 77 —)• 77 process obtained using MC simulation 

is (14.6 ± 0.7)%. The quoted error is systematic, estimated using our work on the 
measurement of the e^e“ —)■ cross section. 

The main background sources for this decay mode are the processes ^ 777 —)■ 37 r °7 
and e’''e“ — )■ Their cross sections were measured in Refs. jlS, 1^, 21 , 22| and in this 

paper (see Sec. m- The number of background events from these sources is calculated to 
be 1.3 ± 0.3 and 0.4 ± 0.1, respectively. It should be noted that the number of data events 
with X^ttOttO ^ equal to 13 is in good agreement with the background prediction based 
on MC simulation: 12 ± 2 for e’''e“ — )■ 777 and 3 ± 1 for 
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D. Decay chain 77' —)• rj Stt^ 

For this decay mode with ten photons in the hnal state there is no background from 
e’''e“ annihilation. The main source of background is cosmic-ray showers. We select events 
containing nine or more photons and no tracks in the drift chamber. The photons must 
have the transverse energy distribution in the calorimeter consistent with the distribution 
for an electromagnetic shower. The parameters F^cai and Pcai must satisfy the conditions 
(|3]). The muon system signal is required to suppress cosmic-ray background. No data events 
are selected after applying these criteria. The detection efficiency for this decay mode is 
( 22.6 ± 1 . 1 )%. 


IV. UPPER LIMIT FOR i e+e” DECAY 


The visible cross section for the process e’''e —)■ rj' is calculated as follows 


N 

exp ’ s 

= 


( 4 ) 


Y.L^ei 

where Ng is the sum of experimental events selected in the hve decay modes, Si is the detec¬ 
tion efficiency in the mode i, which includes the branching fractions for the corresponding r]' 
and 7] decays. The integrated luminosity Lj is equal to Lee = 2.91 pb“^ for the decay modes 
with charged particles and = 2.82 pb“^ for the multiphoton modes. The denominator 
in the formula (jl]) can be represented as LeeSs- For the selection criteria described in the 
previous section Eg = ( 6.2 ±0.4)%. Since the number of selected data events is equal to zero, 
we set the upper limit on the cross section. The technique of Cousins and Highland fol¬ 
lowing the implementation of Barlow 19| is used to calculate the limit with all uncertainties 
included {Ng < 2.32 for 90% CL): 


< 12.7 pb at 90% CL. 


( 5 ) 


The limit on the cross section is translated using Eq. ([2]) to the upper limit on the rj' electronic 
width 

< 0.0020eV at 90% CL. ( 6 ) 

As a test, we perform measurements of the cross sections for the processes 
7i~^7r~7T^, e'^e~ Tr^vr^y and rij ^ The process —)■ 777 was studied in 
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the seven-photon hnal state. The events of these processes are selected with criteria similar 
to those described in the previous section. The obtained Born cross sections a(e~^e~ —)■ 
7r~^7r~7r^) = 11.7 ± 0.2 nb, a(e~^e~ —>■ Tr^Tr^j) = 285 ± 21 pb, a(e~^e~ —)■ 777 ) = 244 ± 30 pb 


are in good agreement with the results of the previous measurements 11.33 ± 0.64 nb 20| 


for e+e —)■ tt+tt 7 r°, 242lg7 pb 21| and OOO^gg^ pb 22| for e+e —)■ 7 r° 7 r° 7 , 300 ± 110 pb 
and 390l}io Q for e+e“ 77 . 
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V. CONCLUSION 


The search for the process e’''e“ —)■ rj' has been performed in the experiment with the 
SND detector at the VEPP-2000 e’''e“-collider. To reconstruct the r/'-meson the hve decay 
chains have been used: rj' —)■ ^vr+vr” followed by the 7 decays to 77 and 37 r°, and 7 ' —)■ 
with Tj decays into 7 r+ 7 r“ 7 r°, 77 , and 37 r°. No data events of the e+e“ —)■ 7 ' process have 
been found. Since the visible cross section for the process under study is proportional to the 
rj' electronic width, we set the upper limit 

P^/^e+e- < 0.0020eV at 90% CL. (7) 

The obtained limit is slightly better than the limit set recently in the CMD-3 experiment 
P^^e+e- < 0.0024 eV je]. 

Using the formula (jl]) we combine the SND (0 events, Lj = 2.91 pb“^, £* = (6.2 ± 0.4)%) 
and CMD-3 (0 events, L* = 2.69 pb“^, Si = (5.3 ± 0.3)%) data and obtain the combined 
upper limits on the electronic width 


T^,^e+e- < 0.0011 eV at 90% CL. 


and the branching fraction [P.,,/ = (0.198 ± 0.009) MeV 10|] 


B{r]' e+e") < 5.6 x 10“^ at 90% CL. 


( 8 ) 

(9) 


The obtained ui; 
predictions 


ed ijp 

aff 


per limit is most stringent but still 30-50 times larger than theoretical 
made in the framework of the Standard Model. 
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